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Introduction. The solution to the problem of transformation of 
phosphorus, nitrogen and silicon forms is studied. This prob- 
lem arises under modeling phytoplankton dynamics in shal- 
low-water bodies including the Azov Sea. The phytoplankton 
dynamics model is formulated as a boundary value problem 
for the system of diffusion-convection-response equations and 
takes into account the absorption and release of nutrients by 
phytoplankton, as well as the transition of nutrients from one 
compound to another. To calculate the initial conditions and 
parameters of the equations under which the steady-state re- 
gime occurs, the software is developed, which is based on the 
model describing changes in phytoplankton concentrations 
without considering current effects. This model is represented 
by a system of inhomogeneous differential equations. Based 
on the developed software, the initial conditions and parame- 
ters of the phytoplankton dynamics model in the Azov Sea are 
calculated experimentally. 

Materials and Methods. A 3D model of phytoplankton dynam- 
ics 1S considered taking into account the transformation of 
phosphorus, nitrogen and silicon compounds based on the 
system of nutrient transport equations. The case of a spatially 
uniform distribution of substances is considered to specify the 
parameters of the model at which the stationary modes occur. 
Because of simplification, a system of ordinary differential 
equations solved through the Runge-Kutta method is obtained. 
Research Results. The software is developed to specify the 
initial conditions and parameters of the phytoplankton dynam- 
ics model considering the transformation of phosphorus, nitro- 
gen and silicon compounds. 

Several numerical experiments are performed under the as- 


sumption that the development of phytoplankton is limited by 
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Beedenue. Padota TocBalleHa HUCCeOBAHHIO pellleHHsA 3ajla- 
4H TpaHchopmamuu dopm docdopa, a30Ta u KpemHua. Jan- 
Had WpoOsiemMa BO3HHKaeT MpH MOJeuHMpoBaHuu jWWHaMUKU 
(UTOMIAHKTOHA B MCJIKOBOJHBIX BOJOCMax, B TOM 4HCIIe B 
A30BCKOM Mope. Mojenb TMHaMuKH (UTONMAHKTOHA cop- 
MyJIMpoBaHa Kak KpaeBad 3ayjada JJId CHCTeMbI ypaBHeHuit 
muddy3uu-KOHBeKUMH-peakIIMH WU YAMTHIBAeT NOTMOMWeHve U 
BbIJ[eICHHe MMTATeCIbHbIX BeEL[CCTB (PUTOMAHKTOHOM, a TakKxKe 
TlepexOJ,| MWHTAaTeCJIbHBIX BeIIeCTB H3 OAHOTO COeMHeEHHA B 
mpyroe. Jia pacueTa HauasIbHbIX yCJIOBHN UM MapaMeTpoB 
ypaBHeHHi, Mp KOTOpbIxX HacTyMaeT CTalMOHapHbIN pexuM, 
pa3paOoTaHo TporpaMMHoe oOeciieyeHve, OCHOBOM KOTOPOTO 
NOcHyKUIAa MOJelIb, OMMChIBAIOMIad H3MCHeCHHA KOHIICHTpa- 
WH MUTOWIAHKTOHa Oe3 yueTa BIIMAHUA TeyeHn. JlaHHas 
MOJeIb IIpecTaBIeHa CHCTEMOM HeEOJHOPOAHBIX OOBIKHOBCH- 
HbBIX TuddepeHiwalbHbIX ypaBHeHui. Ha ocHoBe pa3pavo- 
TaHHOrO TporpaMMHoro oOOecre4eHHA IKCICPHMCHTAJIbHbIM 
00pa30M paccuHTaHbl HayaJIbHbIe YCIOBHA HU MapaMeTpbI MO- 
eu JMHAaMUKH (PUTOMIAHKTOHA B A30BCKOM Mope. 
Mamepuazoi u memoovi. PaccMaTpuBaeTcA TpexMepHad MO- 
eb TMHAaMUKH (UTOMIAHKTOHa C y4YeTOM TpaHcdopmalnu 
coeyMHeHHi docdopa, a30Ta HW KPeMHHA, OCHOBaHHad Ha CH- 
cTeMe ypaBHeHHui TpaHcnopTa OMOreHHbIX BelllecTB. Jia 
YTOUHeCHHA WapaMeTpOB MOJeIu, pw KOTOPbIX HacTyMaroT 
CTal|MOHapHble pexKUMBI, paccMaTpHBaeTcaA Clryyait IpocTpaH- 
CTBCHHO-OJHOpOHOTO paciipeszeneHua cyOcTtanHunui. B pe- 
3yIbTaTe YMpOleHuA MONyyeHa CuCTeMa OObIKHOBCHHBIX 
muddepeninalbHbIx ypaBHeHHH, KOTOpad pellieHa MeTOOM 
Pyure-KyTTbl. 

Pe3yibmamol ucciedoeanus. Pa3zpad0oTaHoO mporpaMMHoe 
oOecnieyeHve Jd YTOUHCHHA HavasIbHbIX yCJIOBMM U Tapa- 
MeTPOB MOJIeJIM JHHAaMUKU ()UTOMIAHKTOHa C yYeTOM TpaHc- 
(bopMallMu coesqMHeEHuM Mocdopa, a30Ta uv KpemHua. IIpose- 
JICHbI HECKOJIBKO YMCJICHHbIX IKCIICPUMCHTOB B TpeTO0%*Ke- 
HHH, YTO pa3sBHTHe MUTOMIAHKTOHa JIMMUTUpyeTcA eMH- 
CTBCHHbIM OHOreHHBIM BelIeCTBOM. B pe3ysIbTaTe BbIYHCJIM- 
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a single biogenic substance. As a result of the computational 
experiment, it can be seen that with the obtained values of the 
initial concentrations and parameters of the equations, station- 
ary modes occur for the system of ordinary differential equa- 
tions describing the case of the spatially uniform distribution 
of substances. 

Discussion and Conclusions. The mathematical model of the 
transformation of phosphorus, nitrogen and silicon forms in 
the problem of phytoplankton dynamics is studied. Stationary 
modes for the system of ordinary differential equations are 
obtained, for which the values of the system parameters and 
initial conditions are determined. The results obtained can be 
used in further simulation of the phytoplankton dynamics con- 
sidering the transformation of phosphorus, nitrogen and silicon 
compounds with account for convection-diffusion, salinity, 


and temperature. 


Keywords: phytoplankton, phosphorus, nitrogen, silicon, bio- 
gen, chemical-biological source, convection-diffusion- 
response equation, Cauchy problem for system of ordinary 
differential equations, stationary mode. 
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TCJIbHOrO IKCIECPHMeHTAa BUHO, 4YTO IPH MOJYICHHBIX 3HAave- 
HHAX HaYaIbHbIX KOHI[CHTpallMu WU MapaMeTpOB ypaBHeHuii 
HacTyliatoT CTallMOHapHble Pe?KHMBI JIJId CHCTCMbI OObIKHO- 
BCHHBIX uddepeHuMaNbHbIX ypaBHeHH, OMMMCbIBAFOIIeH 
ciydai lpOCTPpaHCTBeHHO-paBHOMepHOro _ paciipewesIeHust 
cyOcTaHIUni. 

Oécyacodenue u 3aknrouenus. B padote uccieyqoBaHa MaTeMa- 
THYeCKad MOJIelIb TpaHcdopmMalnu dopm docdopa, a30Ta u 
KPeMHHAl B 3ayjade JMHaMMKH (UTOTAHKTOHA. 

IlonyueHbI cTal[wOHapHble pexKUMBI JJId CHCTeMbI OOBIKHO- 
BeCHHBIX JM@depeHHaIbHbIX yPaBHeHHH, JIA KOTOPbIX ompe- 
JesIeHbI 3HAYeCHHA TapaMeTPOB CHCTeMbI HW HadasIbHble yCJIO- 
Bua. [loryaeHHble pe3yIbTaTbI MOFyT ObITb HCIIOJIb30BaHbI B 
ipollecce MasbHeiiiero MOeIMpoBaHuA JMHaMHKH MuTO- 
IIIaHKTOHa C yueTOM TpaHcdopmallun coeyuHeHUuu docdopa, 
a30Ta HW KPpeMHHA C YYeTOM KOHBEKIMU-7uddy3uu, COeHO- 
CTH, TeMMepaTypbl. 


Ku1r0ueBble C10Ba: (UTOMIAHKTOH, Mocdop, a30T, KPeMHUH, 
OvOreH, XHMUKO-OMOJIOrMYeCKHH HCTOUHHK, ypaBHeHie KOH- 
BeKIMH-Tuddy3uu-peakiun, 3ayaya Komm ia CHCTeMBI 
OObIKHOBeCHHBIX TH@MepeHiwaIbHbIX ypaBHeHH, CTallMoHap- 
HBIM pexKUM. 


Oopazeu Ona uumupoeanua. benosa, tO. B. UccneqoBpanue 
CTaI[MOHapHBIX PelIeHHi 3aaqW TWMHAMUKU (PUTOTJIAHKTOHA C 
yueTOM TpaHcdopmMaliwu coeqwHeHuH dochopa, a30Ta u 
kpemuusa / IO. B. benosa [u gp.] // Becrauk JJonckoro roc. 
TexH. yH-Ta. — 2019. — T.19, Nel. — C. 412. 
https://do1.org/10.23947/1992-5980-2019-19-1-4-12 


Introduction. Because of the development of major cities on the coast of shallow water bodies and river systems 


that flow into these water bodies, eutrophication has become more frequent. The growth of algae in reservoirs is caused 
by an increase in the flow of nitrogen and phosphorus compounds from the adjacent land areas. Each water body is 
unique and requires a thorough study. Field investigations [1] and mathematical modeling are used to explore water 
bodies. Without downplaying the role of field experiments, we can say that mathematical modeling is less costly, and it 
allows us to predict the behavior of the ecosystem. 

To study the Sea of Azov, a three-dimensional model of hydrodynamics [2, 3] including the equations of motion 
in three spatial directions was developed. In [4], this model was made for the case of dynamic rebuilding of the compu- 
tational domain geometry due to the tidal effects. The investigation of this model accuracy is given in [5]. In [6—8], the 
reconstruction of an ecological catastrophe that occurred in 2001 caused by an excessive concentration of algae in the 
eastern part of the Sea of Azov is given. In [9], methods of controlling the suffocation phenomena arising in the Sea of 
Azov were proposed. The [10—12] papers are devoted to studying the dynamics of phyto- and zooplankton. 

The water condition in shallow water bodies is changing rapidly, and mathematical models need to be refined. 
The parameters determination of the three-dimensional model of the phyto- and zooplankton dynamics is laborious; 
therefore, it is proposed to use a simplified model to calculate these parameters. 

The work objective is to improve the parameters of the model of the phytoplankton dynamics considering the 
transformation of phosphorus, nitrogen and silicon compounds, under which stationary regimes occur with the assump- 
tion of a spatially uniform distribution of substances. 

Materials and Methods. The model is based on the system of equations for the transport of nutrients [15, 16], 


the form of which for each F model block is 


Mechanics 


http://vestnik.donstu.ru 


ON 


Vestnik of Don State Technical University. 2019. Vol. 19, no. I, pp. 4-12. ISSN 1992-5980 eISSN 1992-6006 
Becmuuk Jlonckozo 2ocyOapcméennozo mexnuyeckozo ynueepcumema. 2019. T. 19, Noe 1. C. 4-12. ISSN 1992-5980 eISSN 1992-6006 





i aera Oh Oh = div(keradg,)+R., (1) 
Ot ax oy Oz : 


where g,is concentration of the i-th component, [mg/l]; ieM, M={F), FP, F'3, POz POP, DOP, NO3, NO2, NHy, 


Si}; {u,V, w} are components of the velocity vector of the water flow, [m/s]; & is turbulent exchange coefficient, [m/s]; 


R, 1s function-source of nutrients, [mg/I - s]. 
In equation (1), i index indicates the type of substance (Table 1). 
Table 1 
Biogenic substances in model of phytoplankton dynamics 


PO, 
Pop | speded organic phosphor 


por solubie organic phosphor 
NO. 
NO, 
NE, 





Chemical and biological reactions are described by the following equations R, =C,(1—Kyr)dp —Krpdp —KreWs » 


3 
Roop = LSrK 04 — B pn pop NI nop > 
3 
ae — LS pK rer, + K pV pop 7 KS h@nop ’ 
3 
Reo, = LSC, a -1)q, + Koy pop FR pop ’ 


Ga. 3 
Fr (4y04no4on,) +28, (ke + Ky )4 —Kodua, » 


3 
Ruy, = LSC (K,,-1) 
BY (Ivor Inoe2 Ion.) Ivo 


=K, 
f N (Gro, ’ Ao, ’ dv, ) Ano, - no, 


234 no, » 


q, + K, 


2 nit, 


3 
Ryo, = L Snr (K,.. -1) 


(1) 
five = vs.C. (Kee 1) £2 fron) 1 no, 


| 9g +K 
= i N Ce »WIno,» Inu, ) Ino, * Ino, 


23 dno, ’ 


R, =5,C, Ce -1)q, +8 K np , 

Here, K,, is specific breathing rate of phytoplankton; K,,,,is specific die-off rate of phytoplankton; K,,, is 
specific excretion rate of phytoplankton; K,,, is specific rate of POP autolysis; K,,, is coefficient of POP phosphatifica- 
tion; K,,,, is coefficient of DOP phosphatification; K,, is specific rate of ammonium oxidation to nitrites under nitrifica- 
tion; K,,1s specific rate of nitrite oxidation to nitrates under nitrification; s,, s,, s,,are normalization coefficients be- 


tween the content of N, P, Si in organic matter. 
Phytoplankton growth rate is determined by the following expressions: 


Ci. = Kup, min fe (Gro, ) a Cee Inn, )} 


C;, Sg min fe Ce ee (edatde, ne (ds ) 4 


where Ky is maximum specific phytoplankton growth rate. 
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Fig. 1. Model scheme of biogeochemical transformation of phosphorus, nitrogen and silica forms 
Functions describing nutrient content 


1 po, 


— for phosphorus: f,(g,,. )= 
(dr Ino, + Ko, 


where K,,, 1s phosphates half saturation constant; 
1 Is; 
— for silica: f, (45,)=—_ 
9a t Ks 


where K,. is silica half saturation constant; 


— for nitrogen: f, (Gro, Teste ) = fe Ce »WIno,» Inn, ) +90) (Cae ) ) 


(1) (Ino, ale Ayo, ) exp(—K dT.) (2) 
tn CR Reet ee See Ix CE =e 
Kyo, + (no, + Ino, ) Kon Cap, 
where K,,, is nitrates half saturation constant, A,,, 1s ammonium half saturation constant, A,,, 1s ammonium inhibi- 
tion ratio. 
For the system (1), it is necessary to specify the vector field of the water flow velocities, as well as gq. initial 


values of the concentration functions 
q, (x, y,z,0) =G (x,y,z), (x,y,z) eG , t=0,ieM. (2) 
Assume % boundary of G cylindrical region is piecewise smooth and XY =X, UX, Uo , where &,, is surface of 
the reservoir bottom, X, is still water surface, o is lateral (cylindrical) surface. Suppose uw, is vector component of the 
water flow velocity normal to 2 , and 7 is outward normal vector to X . Forg, concentrations, we assume: 


— on side boundary: 
g,=9, ono,if u, <0,ieM; (3) 


Oq. 
AG, og, 


=0, on o,if u, 2>0,ieEM; 
On On 


—on X, there is reservoir surface: 


OTe G. gen (5) 


IZ 


—on the bottom %,, : 


Og. ; 
ke = 6,,9¢,. t€{Fi, Fo, Fs}, 
Zz 
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k at =«,,g,, i€ {PO4, POP, DOP, NO3, NO>, NHy, Si}, (6) 
IZ 


where &, ,, €,,are sedimentation rates of algae and nutrients to the bottom. 





Stationary Mode. Consider the case of a spatially uniform distribution of substances (phytoplankton, forms of 
phosphorus, nitrogen and silica); then each of the equations (1) is simplified; and as a result, we get the following sys- 


tem of ordinary differential equations (ODE): 


dq, a 
dt = C, (l— Keds a Kode — Kaede 0H 195 (7) 





oe 
—_ => s pK ED qr, = Ky drop — By pop» 


At 
—_— => S»K BE q a K, pd pop ~ or pop 9 


7 
——. m= 9's,C (K, 2a, +K, vA pop + K, 1 pop » 

Ayn, -¥s.C =) fy (NH,) Ys, (Key + Kx )4y —Kud 
i a OO ENO .NO.NH) " a 2 Ue eee 
dq 3 (NO,, NO,, NH) q 
Se Di SG. (Kaeo a 1) fx (NO,NO,»NH,) Gr g + Ky yn, —K3.9yo0,> 
dt i=l I (NO,, NO, , NH,) Ano, T Ino, 

d ") (NO,, NO,, NH 

do, =>s,C, (Ky, 1) NO NO NH) doy qd, EK Gn : 

dt i=l fy (NO,, NO,,.NH, ) Avo, * Ino, 





“ = 5,C, Ca -1) 4, gore ae ae 


We solve the system of ordinary differential equations by the Runge — Kutta method [15-17]. We will conduct 
several numerical experiments, assuming that the development of phytoplankton depends on a single limiting substance. 


Research Results. For the ODE system (7), we calculate the initial conditions and parameters of the equations 


at which the stationary regimes occur. Let us take the initial concentration values: g, (0)=2.5 mg/l, q,, (0) = 2.6 mg/, 
d,,(0)=0.91 mg/l, dpop(0)= 9.07 mg/l, yop (0)= 9.07 mg/l, Apo, (0) = 0.005 mg/l, qy,,(0)= 0.11 mg/l, yo, (0) 
= 0.0178 mg/l, qyo, (0) = 0.304 mg/l, g,,(0)= 0.4 mg/l; coefficients: K,,. = 2.8 day", K,,=0.15 day”, K,,,= 0.05 day” 
. Kg = 0.15 day" , Kpp = 0.015 day" , Kpy = 0.02 day" , Kpn = 0.1 day" , Kap = 0.9 day” , Ko3 = 2.5 day“, Kysi = 1.46 


day, sp=0.01, sy = 0.016, ss; = 0.023, Kyo, = 0.024, Ky, = 3.0, Ky, =2.0, K, = 3.0. 


The obtained stationary modes of the ODE system (7) on the assumption that the development of phytoplank- 
ton is limited by a single nutrient (phosphorus, nitrogen or silica) are shown in Fig. 2-4, respectively. Fig. 2 describes 
the effect of phosphorus on the development of various phytoplankton species; Fig. 3 describes the effect of nitrogen on 


the development of various phytoplankton species; Fig. 4 describes the effect of nitrogen on diatom development. 
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chy, ms 
2.6 

















Fig. 2. Stationary mode of ODE system under assumption that phytoplankton development is limited by phosphorus: 
a) green algae (ChV), b) green-blue algae (AF-A), c) diatom (SC), d) suspended organic phosphorus (POP), e) soluble 
organic phosphorus (DOP), f) phosphates (PO,) 
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Fig. 3. Stationary mode of ODE system under assumption that phytoplankton development is limited by nitrogen: 
a) green algae (ChV), b) green-blue algae (AF-A), c) diatom (SC), d) ammonium (NH,), e) nitrites (VO>), f) nitrates (VO3). 


sc, mg/l 
10 


Fig.4. Stationary mode of ODE system under assumption that phytoplankton (diatoms) development is limited by silica: 
a) diatom (SC), b) silica (Si). 
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The result of the computational experiment shows that with the above values of the initial concentrations and 
parameters of the equations, stationary modes occur for the ODE system (7), which describes the case of a spatially 
uniform distribution of substances. The obtained values will be used in further simulation of the spatially inhomogene- 
ous distribution of substances, saltiness and temperature considering the movement of the aquatic environment [18]. 

Conclusion. A mathematical model of the transformation of forms of phosphorus, nitrogen and silica in the 
problem of phytoplankton dynamics is studied in the paper. The case of spatially uniform distribution of substances 
(phytoplankton, forms of phosphorus, nitrogen and silica) is considered. The system is divided into three systems of 
ordinary differential equations, each of which simulates the dependence of phytoplankton growth on a single nutrient. 
These systems are solved by the Runge-Kutta method (Fig. 2—4); stationary modes are obtained, for which the values of 
the system parameters and initial conditions are determined. 

The results obtained will be used for the further simulation of the phytoplankton dynamics considering the 
transformation of phosphorus, nitrogen and silica compounds, taking into account diffusion-convection, saltiness, and 
temperature. 
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